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We study Holstein polarons in three-dimensional anisotropic materials. Using a variational exact 
diagonalization technique we provide highly accurate results for the polaron mass and polaron radius. 
With these data we discuss the differences between polaron formation in dimension one and three, 
and at small and large phonon frequency. Varying the anisotropy we demonstrate how a polaron 
evolves from a one-dimensional to a three-dimensional quasiparticle. We thereby resolve the issue of 
polaron stability in quasi-one-dimensional substances and clarify to what extent such polarons can 
be described as one-dimensional objects. We finally show that even the local Holstein interaction 
leads to an enhancement of anisotropy in charge carrier motion. 

PACS numbers: 71.38.-k,71.38.Cn,71.38.Ht 



I. INTRODUCTION 



In materials with a local coupling of the charge carriers 
to optical phonons a tendency towards formation of po- 
larons prevails. How polarons form, and which properties 
they possess, has been intensively studied using the Hol- 
stein model (for a recent review see Ref . [l]) . Much work 
was devoted to the question how polaronic properties de- 
pend on dimensionality. For zero phonon frequency po- 
laron formation is known to be fundamentally different 
in dimension one (ID) and three (3D): While in ID a po- 
laron exists for all coupling strengths, in 3D a transition 
from a free electron to a self-trapped polaron occurs at a 
finite coupling strength^. While the 3D polaron is small, 
in ID large polarons with a size of many lattice constants 
may exist 3 . This qualitative difference is peculiar to the 
static case. At finite phonon frequencies no strict self- 
trapping transition occurs^. Instead a smooth crossover 
from a light electron to a heavy small polaron takes place 
with increasing coupling strength. Only in the adiabatic 
regime of small phonon frequencies, large and heavy po- 
larons can be found in 1D S . In the anti-adiabatic limit 
of large phonon frequencies polaron properties do not 
depend on dimension. For intermediate phonon frequen- 
cies numerical studies^ reveal that in 3D polarons form as 
smaller and heavier quasiparticles if compared to the ID 
case, but a clear distinction as for zero phonon frequency 
or in the adiabatic regime is not found. 

Most work on the Holstein polaron problem was per- 
formed for ID or isotropic systems. Much less atten- 
tion is paid to the question how polaron properties in 
anisotropic materials, e.g. molecular crystals 7 , interpo- 
late between the ID and 3D behaviour. In particular, 
even polarons in quasi- ID systems, such as in organic 
conductors like the Bechgaard salts, are frequently de- 
scribed by the ID Holstein model, neglecting the possibil- 
ity of electron motion transverse to a quasi- ID chain. For 
zero phonon frequency Emin showed^ that already tiny 
electronic transfer integrals perpendicular to the chain 



suffice to destabilize a large ID polaron, resulting in a 
free electron as in 3D below the self-trapping transition. 
This result not only puts the existence of large (adia- 
batic) polarons in real materials in doubt, but also raises 
the question to which degree polarons in quasi-lD sys- 
tems can be described as true ID polarons. In the present 
work we study to which extent the scenario described for 
zero phonon frequency translates to finite phonon fre- 
quencies beyond the extreme adiabatic regime. 



II. ANISOTROPIC POLARON MODEL 

To address this fundamental question, we consider the 
3D Holstein polaron model with anisotropic hopping 



H = - E ^E^+e, 

i—x,y,z n 



E^n + M c n c n + ^0 E b lX 



(1) 



Here n labels the sites of a cubic lattice, and e$ for 
i = x,y,z denotes the unit lattice vector in the re- 
spective direction, ojq is the phonon frequency, and e p 
the electron-phonon coupling strength. For the elec- 
tron transfer integrals, we set t x = iii, t y = t z = t±. 
As t±/tu grows, the system described by the Hamilto- 
nian |T]) evolves from a ID (t±/t\\ = 0) to an isotropic 
3D system (t±/tu — 1). For intermediate values of t±/tu, 
the system is anisotropic but symmetric with respect to 
the y, z-direction, i.e., the chains in a quasi-lD system 
(tx/~t\\ <C 1) are oriented along the x-direction. 

In addition to Ref. |U, the anisotropic Holstein 
model (Q} has been studied in Refs. IMlol in certain 
limiting cases, or using approximate variational tech- 
niques. At present, no conclusive answers for small-to- 
intermediate phonon frequency have been given. In this 
work, we present accurate results for the infinite system 
using the variational exact diagonalization technique de- 
veloped in Ref. [ll|. The error of all results is generally 
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FIG. 1: (Color online) Polaron radius R as a function of cou- 
pling A, for dimension D = 1 (t± = 0) and D = 3 (isotropic 
case tu = tj_). The dotted curve gives the radius for lu — 
and ID, which is defined according to Eq. ([2]) using the elec- 
tron density ^(r)! 2 instead of the correlation function x( r )- 



much less than 1%, and about 1% for the polaron radius 
in the crossover regime at small luq in 3D. Our presen- 
tation concentrates on small-to-intermediate phonon fre- 
quencies wo A|| = 0.1 . . • 2.0. We justify this choice in the 
next paragraphs. 



III. NUMERICAL RESULTS AND DISCUSSION 
A. ID and 3D cases 

The destabilization of a polaron at small £j_A|| oc ~ 
curs only if the polaron is large.® We therefore begin 
our discussion with the basic question in what cases 
large polarons exist at finite phonon frequencies. Us- 
ing the electron-phonon correlation function x(r) = 

Env^olOi+r +&n+r)cj,c n |'0o)i the polaron radius, which 
measures the extension of the phonon cloud surrounding 
the electron, is defined as 



£ r (^) 2 x(r) 



-,1/2 



2£rX(r) 



(2) 



and similarly for R± with r x replaced by r y or r z . Note 
that for the Holstein model, £ r x( r ) = 2y / £^7^o- 

How the distinctive features of polaron formation 
found for zero phonon frequency persist at small WoA|| i s 
manifest in a comparison of the ID (£j_A|| = 0) and 3D 
(£_i_A|| = 1) curves for R = R\\ = R± in Fig. [TJ given as 
a function of A = e p /2Dtn. The steep decrease of R in 
the vicinity of A = 1 for 3D, uio/t = 0.5 is a precursor 
of the self-trapping transition for loq — 0, which takes 
place at AjP w 0.90. In ID the polaron radius decreases 
steadily with A, and no crossover region can be identi- 
fied. Although the latter behaviour is reminiscent of the 
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FIG. 2: (Color online) Polaron mass m as a function of cou- 
pling A, for the parameters as in the previous figure. 



ID large-to-small polaron crossover at luq = 0, we ob- 
serve that even for small phonon frequency woA|| = 0.1 
no large polaron is found. In Ref. [l_3j it was shown that, 
independent of dimension, for weak coupling the polaron 
radius approaches 



(3) 



(compare the ID and 3D curve for WqA|| = 0.5 in Fig.[T]). 
The polaron radius is bounded by that limiting value, 
e.g., for wo A || = at mos t R\\ < 3.2. To obtain a 
large polaron of radius greater than 5 or 10, wo/£|| must 
be less than 0.01 to 0.04 for weak coupling e p /wo, A <C 1. 
Furthermore the polaron shrinks rapidly with increasing 
coupling. Even for ujq = a coupling of, say, A > 0.25 
leads to a polaron radius of only R < 2. The condition on 
the phonon frequency this becomes even more restrictive 
if we ask for large polarons with substantial mass, in 
contrast to the light large polarons at weak coupling. 

The polaron mass, which is probably the most impor- 
tant quantity to characterize a polaron, is defined as 



d 2 E{k) 



dkf 



k=0 



(4) 



where E(k) is the groundstate energy at momentum k. 
We set mil = m„ and m± — m y — m z . In the absence of 
electron-phonon coupling (e p = 0), the masses are given 



by m{ = l/2i|| 



= 1/2£_l. At lu = the (3D) 



transition from a free electron to a self-trapped polaron 



coincides with a jump of m; from rrii 



(below the 



transition) to m, = oo (above the transition). For finite 
wo, m depends analytically on the coupling strength e^. 

In Fig. [5]we show the polaron mass m — mu as a func- 
tion of the coupling strength A = e p /2Dt, for different 
luq in dimension D = 1,3. Remember that the data for 
m correspond to the exact numerical solution of the Hol- 
stein model. As for the radius the distinctive differences 
of polaron formation in ID and 3D are manifest in a 
comparison of the corresponding curves for wq/£|| = 0.5. 
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While in ID the polaron mass increases steadily with A, 
a steep increase of my occurs in the vicinity of for 
3D. One can therefore locate, for small phonon frequency 
u>o/tu < 1 in 3D, the crossover from a light particle to 
a heavy polaron in the region A ~ 1, while in ID no 
crossover region can be identified. The differences be- 
tween ID and 3D vanish for larger phonon frequencies 
LUo/t\\ ^ 1 (uot shown here), where the mass indicates the 
continuous evolution from a light to a heavy polaron with 
increasing A. In the anti-adiabatic regime u>o/tu 3* 1, the 
relation mi/m® = exp(<? 2 ) with g 2 — s p /loq holds inde- 
pendent of dimension. 

The behaviour of the ID polaron mass for small 
phonon frequency is reminiscent of the large-to-small po- 
laron crossover in the coo = 0- limit. However, for ujq = 
the mass of the ID polaron is strictly infinite for any cou- 
pling strength in ID. The ujq — 0-result therefore does not 
provide a quantitative prediction for the polaron mass at 
any finite phonon frequency. Owing to that restriction, 
it is also not sufficient to establish the existence of large 
polarons with substantial mass renormalization even for 
the adiabatic case Wo/tii <C 1. Instead we see in Fig. [3J 
that even for a small phonon frequency woA|| =0.1 the 
mass of the ID polaron is still close to unity for those 
coupling strengths A < 0.25, for which a large polaron 
is found for sufficiently small u>o/t\\. To obtain polarons 
with large mass and large radius therefore requires small 
coupling strengths A < 1 (to allow for a large radius) 
and even smaller phonon frequencies wo/^|| ^ A (to al- 
low for a large mass) . Consequently, such heavy adiabatic 
large polarons may exist only at extremely small phonon 
frequencies Wo/iii < 0.01. We note that (i) in general 
no large-to-small polaron crossover occurs even at small 
phonon frequencies in ID, (ii) the adiabatic regime of 
large (heavy) polarons occupies only the very narrow pa- 
rameter region of tiny phonon frequencies, thus (iii) the 
destabilization of large quasi-lD polarons at small t±/t^ 
Emin obtained for ujq = can be expected to be rele- 
vant only in that narrow region. Qualitatively, we ex- 
pect strong deviations from the findings for luq — even 
for small-to-intermediate phonon frequencies beyond the 
extreme adiabatic regime. 



B. Spatial anisotropic case 

Let us now begin with the quantitative discussion of 
the anisotropic case. Since in the extreme adiabatic 
regime of large (heavy) polarons the question of polaron 
properties in anisotropic materials was already settled 
by Emin£ we exclude this narrow region from our further 
consideration. Instead we concentrate on the broader 
regime of small-to- intermediate phonon frequencies. 

Obviously, the polaron radius does not indicate per 
se to which degree the polaron ceases to be a true ID 
particle for small tx/tu. We may instead use the ratio 
Rx/R\\ (see Fig. [3]), which we interpret as a measure of 
the "one-dimensionality" of the polaron: For Rx/R\\ = 



t 1 1 1 1 1 1 r 




FIG. 3: (Color online) Ratio Rx/R\\ of polaron radii as a 
function of anisotropy tx/tu, for uo/t^ and e p /t^ as indi- 
cated. Here the dotted curve gives the weak-coupling result 
Rx/R\\ = y/tx/t\\- 



(Rx/R\\ = 1), the polaron is a fully ID (isotropic 3D) 
object. We see that already for small tj_A|| the ratio 
Rx/R\\ deviates significantly from zero, but Rx/R\\ is a 
continuous function of tx/tu in contrast to the behaviour 
in the wo = 0- limit. For weak coupling, we recover 
Rx/R\\ = y/tx/t\\, as can be derived by perturbation 
theory^. For strong coupling, the asymptotic behaviour 
R±/R\\ = tx/tu is approached. For intermediate cou- 
pling (e p /t\\ = 4.5 in the figure) we observe how the 
curve interpolates between the strong-coupling straight 
line at smaller ij_A|| and the weak-coupling square root 
at larger tx/tu- To understand this observation, we shall 
now use the polaron mass to further clarify how the po- 
laron properties change with tx/t\\. 

In Fig. 2] we show the polaron mass my as a function of 
the coupling strength A = e p /2(t\\ +2ij_), for different ojq 
and anisotropy tx/tu. We see that the curves for inter- 
mediate values of tx/tu smoothly interpolate between the 
ID and 3D case. In contrast to the ojq = 0-limit a sud- 
den transition from ID to 3D behaviour at small tx/tu 
is missing. Nevertheless, for woA|| = 0.5 and A > 1, the 
mass my increases by a huge factor if tx/tu grows from 
zero to one. Note that this scenario is opposite to po- 
laron destabilization in the loq = 0-limit, where a large 
polaron evolves into a free electron at weaker coupling 
A < 1 (cf. R{uj = 0) in Fig. [J). 

In Fig. [5] we show toii as a function of tx/tu for fixed 
E p . Generally, the polaron mass my decreases with in- 
creasing tx- For larger phonon frequencies, mu depends 
only weakly on tx/t\\, in agreement with our considera- 
tions above. For small phonon frequency (wo/iy = 0.5), 
we see how a heavy ID polaron with m,||/mjj 3> 1 evolves 
into a light 3D particle with my /mjj ps 1 as tx/t\\ grows 
from zero to one. This behaviour occurs if the coupling 
is sufficiently large to create a heavy polaron in ID (i.e. 
e p /2i|| ^> 1) but below the heavy polaron crossover re- 
gion A ~ 1 in 3D seen in the left panel in Fig. 0] (i.e. 
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FIG. 4: (Color online) Polaron mass m\\ for different 
anisotropy t±/tn, and phonon frequencies u>o as indicated. 
In each panel, eleven curves for tx/t\\ = 0.0 . . . 1.0 in steps 
of 0.1 are shown. Note that for fixed A = £ p /2(tii + 2tj_) the 
coupling e p depends on t±. 
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FIG. 5: (Color online) Polaron mass mu as a function of 
anisotropy tx/tn, for uo/i|| and £ P /tii as indicated. 



e p /6i|| < 1). Exactly then the ratio R±/R\\ interpo- 
lates between the strong- and weak-coupling curve (cf. 
Fig. [3]) . Based on this condition a rough estimate for the 
t±/t\\ value when the particle finally becomes light (i.e. 
rn||/m | ps 1) is provided by the comparison of coupling 
energy to kinetic energy, yielding e p /2(t\\ + 2t±) > 1 or 

*i/*u ~ ( £ p/h - 2 V 4 > e -§- *-l/*ii ~ °- 5 for £ p/h = 4 - 5 - 

This estimate docs essentially not depend on phonon fre- 
quency. Especially, our findings persist also for wo/i|| < 
0.5. 

The observed evolution from a heavy to a light particle 
at small-to-intermediate phonon frequency is the equiv- 
alent of polaron destabilization in the luq = 0-limit. We 
however note significant differences to the ujq = see- 
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FIG. 6: (Color online) my, m± and mass renormalization 
anisotropy A as a function of tx/tu, for £p/t\\ = 4.0. Left 
panel: For t±/t\\ — > 0, the renormalization mx/mx is given 
by the inverse quasiparticle weight 1/Z (dotted line). Right 
panel: In addition, A is shown for a 2D system. 



nario. First, it requires coupling strengths for which the 
ID polaron is not large. If the coupling is weak enough 
to allow for a large ID polaron, no significant change 
of the polaron mass can occur on increasing tx/tn since 
the ID mass is already of the order m||/mjj ~ 1. Recall 
that we exclude the extreme adiabatic regime, where this 
statement may be violated for extremely small phonon 
frequencies, from our considerations. Second, in contrast 
to the behaviour at zero phonon frequency, the change 
in mil takes place in a large range of tx/t\\, and does not 
occur as a rapid crossover in the vicinity of a small tx/tii- 
value. These differences clearly show that the concept of 
polaron destabilization in quasi- ID systems is only of re- 
stricted relevance at finite phonon frequencies. This is a 
consequence of the fact that the arguments for uj = 
rely on the existence of large polarons. As discussed be- 
fore, those fail to exist even at small phonon frequencies 
wo/*|| £0.1. 

A repeated question raised for quasi-lD systems is 
whether a local or isotropic electron-phonon interaction 
can enhance the anisotropy in the electronic transfer inte- 
grals, leading to mx/m\\ > (tx/t^)^ 1 . Large anisotropics 
in polaron mobility that exceed the estimates from, e.g., 
band structure calculations by several orders of magni- 
tude can be explained if acoustic phonons are taken into 
account^?. It is common belief that a similar mechanism 
is absent in the Holstein model due to the purely local 
interaction with isotropic optical phonons. 

As a measure for the anisotropy of mass renormaliza- 
tion, we define 
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■m±/m j_ — m||/mjj 
m»/rm. 



(5) 



or equivalently m±/m\\ = (1 + A)mj_/mjj. In the anti- 
adiabatic strong-coupling limit, with mi/m® = exp(g 2 ), 
the mass renormalization is isotropic (A = 0), but for 
smaller phonon frequency A is possible. An en- 
hancement of anisotropy corresponds to A > 0. 

In Fig. [5] we show mil, mj_ and the ratio A for three 
choices of parameters. We see that generally m± > mu, 
and A is maximal for t± = 0, while of course A = 
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FIG. 7: (Color online) Fractional difference 5 between m\\/m® 
and Z as a function of wo/tii, in ID with t± = 0, for fixed 
g 2 = £ p /uj(), or fixed e^/ty in the inset. 



materials, including the ID and isotropic 3D case. In the 
extreme adiabatic regime, where large heavy ID polarons 
exist for tiny phonon frequencies wo/iy <C 0.1, the sce- 
nario of polaron destabilization given by Emirs— applies. 
We present complementary results for the broader regime 
of small-to-intermediate phonon frequency and small-to- 
strong coupling strength. In combination with Emm's 
work at least the most fundamental answers to the Hol- 
stein polaron anisotropy problem are thereby given. The 
regime of very strong coupling or large phonon frequency, 
where a polaron is always small or its properties do not 
depend on dimension, deserves no further exploration in 
that context. 



for t_i_/i|i = 1. Evidently, A is largest for small phonon 
frequency and large coupling, when deviations from the 
anti- adiabatic limit with A = are most significant. 

From first order perturbation theory in t± we conclude 
that A — ► 6 for t± — > 0, where 



(mS/m||) ti= o - Z 



Z 



(6) 



is the fractional difference between the quasiparticle 
weight Z = |(V'o|cj c=0 |vac)| 2 at momentum k = (\i/jo) 
denotes the k = groundstate) and the inverse mass 
renormalization (mjj/m||) ti= o at t± = 0. This relation 
holds independent of dimension. As a consequence of 
the momentum dependence of the polaron self-energy, 
S ^ is possible 6 . That S > 0, and thus A > 0, re- 
sults from effective long-range hopping processes induced 
by electron-phonon-coupling, which reduce (m®/m\\)t ± =o 
compared to l/Z. In a certain sense, 6 measures the 
deviation from the anti-adiabatic strong-coupling result 
mi/m® — Z^ 1 = exp(g 2 ), when 5 = 0. 

In Fig. [7] we show 6 for fixed g 2 as a function of u>o- For 
u! /t ^S> 1, the ratio S is close to zero as the anti-adiabatic 
limit is approached. For luq — > and fixed g 2 , the ratio 
S converges to zero since the weak coupling regime is 
reached. Decreasing ujq or increasing e p , the ratio 5, and 
thus A, can be made very large. For the parameters plot- 
ted, A can attain values as large as 27%. Note that A is 
already large for small coupling eo/tu provided Wo/tu is 
small (see inset in Fig. [Tj), which implies that the mech- 
anism described here is relevant for realistic materials. 



Our findings have two major implications. First, po- 
laron destabilization occurs only in the extreme adiabatic 
regime, but the concept does not carry over to finite 
phonon frequencies without significant modifications. Al- 
though it is possible to drive a heavy ID polaron into 
a light 3D particle by increasing the perpendicular elec- 
tronic transfer integral, a smooth crossover takes place in- 
stead of an instability. In any case, polarons in quasi-lD 
materials can be described by the ID Holstein model for 
small perpendicular transfer integrals. The second impli- 
cation is that even a local Holstein interaction enhances 
the anisotropy in electron motion. For small phonon fre- 
quencies this enhancement is of significant size already at 
weak-to-moderate coupling. As a question for future re- 
search we mention that the anisotropy enhancement be- 
comes even more important for longer-ranged electron- 
phonon coupling. We stress that the frequent assump- 
tion that the Holstein interaction leads only to isotropic 
changes of material properties is not true. Experimen- 
tally observed anisotropic behaviour provides no indica- 
tion against a dominating short-range Holstein-like cou- 
pling. 
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